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i.  ^Introduction. 

Scattering  of  light  by  molecules  is  commonly  dealt  with  through 
Rayleigh-Debye  theory  Such  approach  is  known  to  be  applicable  to 

molecules  whose  effective  index  of  refraction  is  close  to  unity  and  to 
imply  anuiaber  of  approximations  which  may  have  rather  severe 
effects.  **  In  this  paper  we  face  the  problem  by  adapting  the  nethou 
devised  by  Johnson^*  and  successfully  used  to  calculate  the  electronic 
states  of  large  molecules /^to  the  scattering  of  electromagnetic  waves. 
Accordingly,  we  model  a  molecule  as  a  cluster  of  spherical  scatterers, 
possibly  of  different  radii  and  refractive  indexes.  A  plane 
electromagnetic  wave,  incident  to  the  cluster,  undergoes  multiple 
scatterings  which  we  account  for  by  expanding  the  scattered  wave  in  a 
nulticentered  series  of  vector  harmonics.  The  expansion  coefficients 
turn  out  to  be  the  solutions  of  the  system  of  linear  equations  obtained 
by  expanding  both  the  incident  field  and  that  within  the  spheres  in  a 
series  of  vector  harmonics,  and  imposing  the  boundary  conditions  across 
the  surface.  Due  to  the  presence  of  the  incident  plane  wave,  the  above 
system  is  nonhomo geneous . 

The  method  outlined  above  does  not  require  any  approximation, 
except  for  the  truncations  of  the  expansions  in  vector  harmonics,  in 
order  to  get  a  system  of  finite  order.  For  large  clusters  and  for 
convergency  reasons,  however,  one  might  have  to  solve  rsther  large 
systems.  Anyway,  if  the  cluster  possesses  symmetry  properties,  as  is  the 
case  for  actual  molecules,  one  can  use  group  theoretical  techniques  to 
get  the  system  in  factorised  form.  By  expanding  both  the  incident  and 
the  scattered  field  in  symmetry-adapted  combinations  of  vector  harmonics, 
indeed,  the  system  factorizes  in  ouch  the  same  way  as  quantum  mechanical 
secular  equations,  except  that  the  inhomogeneous  terms  are  not 
independent  of  the  row  index  of  the  multidimensional  irreducible 
representations.  As  a  consequence,  we  have  to  solve  all'  the  systems 
arising  out  of  the  factorization  procedure.  It  will  become  apparent. 
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however,  Chat  this,  peculiarity  of  the  present  problem  does  not  cancel  the 
undoubted  advantages  of  the  group  theoretical  techniques. 

2.  MULTIPOLAR  EXPANSION  OF  THE  FIELD. 

The  cluster  whose  scattering  properties  we  want  to  study  is 
composed  of  N  non-magnetic  spheres  whose  ceuters  lie  at  R  and  whose 
radii  and  (possibly  complex)  refractve  indices  are  ab  and  n  , 
respectively.  We  refer  the  cluster  to  a  fixed  set  of  axes  and  choose  £he 
direction  of  incidence  of  the  incoming  plane  wave  through  direction 
cosines  of  its  wavevector.  A  straightforward  calculation  along  the  lines 
sketched  by  Jackson10  allows  us  to  write  the  field  of  a  circularly 
polarised  plane  wave  ofwavevector  k  a  s  _  _ 


*  £  WnLM(k)[jL(kr)xLM(r)  +  n|vxjL(kr)XU((r)]  (2-1a) 


IB*1*  «  nE^5  (2— 1b) 

-n  -n 

with  n  -  +  1  according  to  the  polarisation  and 

WnUM(k)‘=  +  1ne2)  •  X*M(k)  (2-2) 

e  * 

where  e  and  e  ^  are  unit  vectors  orthogonal  to  k  and  to  each 

other.  ~In  the  above  equations  the  vector  spherical  harmonics,  X  _  , 
are  defined  according  to  Jackson.1 

The  field  scattered  by  the  cluster  is  expanded  in  a  multicentered  series 
of  multipoles  including  only  outgoing  spherical  waves  at  infinity: 


(2-3a) 


(2-3b) 


with  r  -  r  -  R  .  The  superscript  1)  on  the  spherical  Uankel 
functions'*  of  ~the  first  kind  will  be  omitted  throughout  for  simplicity. 
As  regards  the  field  within  the  spheres  we  remark  that  the  present  theory 
is  not  restricted  to  homogeneous  spheres,  provided  the  n  's  are 
spherically  symmetric.  Therefore,  within  the  a  -  th  sphere  we  can 
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where  if  and  s£  are  Che  solutions,  regular  at  r^  «  0,  of  the 
equations1* 


and 


respectively. 
na*s,  Rj(ra) 


L  ®  aaao  a  J 

Of  course  for  unifora _ 

S “  (r  )  5  j  (Kr)  with  K  -  kn  . 
L  a  L  a  a  a  a 


(2-5«) 


(2— 5b) 


3.  SYMMETRY  ADAPTED  MJLTIPOLAR  EXPANSION  OF  THE  FIELD 

The  expansion  coefficients  in  equations  (2-3)  and  (2-4)  are 
solutions  of  the  system  of  linear  nonhonogeneous  equations  obtained  by 
imposing  on  E  and  B  the  customary  boundary  conditions  at  the  surface  of 
each  one  of  "the  spheres.  However,  by  first  exploiting  the  symmetry 
properties  of  the  cluster  and  thus  putting  equations  (2-1)  and  (2-3)  in 
symmetry-adapted  form,  we  shall  get  the  above  system  in  factorized  form. 

First  of  all  we  partition  the  cluster  into  sets  of  spheres  which 
are  transformed  into  each  other  by  the  symmetry  operations.  This  does 
not  imply  renumbering  of  the  spheres  but  only  association  of  each  site 
index,  a,  to  the  appropriate  set  index,  o.  Then  we  construct 
symmetrized  combinations  out  of  the  vector  harmonics  centered  at  the 
sites  of  each  set.  To  this  end  we  remark,  that  if  fL  (kr)  is  a  spherical 
Bessel  or  Hankel  function  then 

M^Xjr)  -  -fL(k?)ll(r) 


Vx/L(kr)XLM(r) 
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the  T's  being  irreducible  spherical  tensors.12  Therefore,  both  the 
magnetic  and  the  .electric  2r  -pole  fields  traishini  according  to  the 
representation  DlL'  of  the  full  rotation  group.12 ,1J*27  Accordingly,  let 
S  be  a  syanstry  operation  such  that  ____ 


5, 


with  a  and  0  ±n  the  sane  set,  of  course,  and  let  0g  be  the  associated 
operator,  then  14  if  S  is  a  proper  rotation 


<Wkra>?u.<r.>  •4^£D^l(S)Xur  (f8) 


(3-2s) 


Vx/JkrjX^r.)  •  £ ’*/t(k',e)o2i(S)Xur (f8). 


(J-2b) 


However,  if  S  is  an  improper  rotation  one  oust  take  into  account  that  the 
parity  of  is  (-)L  while  that  of  7  *  ft  ?  M  is  (-)L**  .  Therefore 

for  improper  rotations,  the  right-hand  side  of  equation  of  (3-2a)  must  be 
multiplied  by  (-)L  and  the  right-hand  side  of  equation  (3-2b) 
by  (-) L*'  and  the  argument  in  o1^  must  be  understood  as  the  proper 
rotational  part  of  S.  As  a  consequence,  if  the  symmetry  group  of  the 
cluster  also  includes  harmonics  belonging  to  the  rows  of  the  v  -th 
irreducible  representation,  one  has  to  apply  the  projection  operators15,16 

•tSOs.o, 

both  to  the  magnetic  and  the  electric  2^fl>oles.  Thus  we  can  write 


(3-4.) 


Kvpo  -  E  E  bvp®  vxh  (krjX  {rj 
'"L  aeofiTe*  L  °“LM  -a 


(3-4b) 


for  the  combinations  of  magnetic  and  electric  2h.  poles,  respectively, 
centered  at  the  sites  of  the  o-th  set.  The  superscripts  v,  p 
indicate  that  the  combination  belongs  to  the  p-th  row  of  the  v-th 
irreduciole  representation  and  the  index  N  recalls,  when  appropriate, 
that  one  can  get  more  than  one  set  of  basis  functions  for  a  given  L.  The 
scattered  field  can  thus  be  written  in  symmetrized  form  as 


£  £[£«  sr  ♦  £  »ss.  c] 


(3-5.) 


the  corresponding  expression  for  iB  8'  being  obtained  through  the 
iiaxuell  equation  iB  ■  -  7xE;  the  inde^  n  denotes  the  polarization  as 
in  the  preceding  section. 

Comparison  of  equation  (3-5)  with  equation  (2-3)  shows  that 

T  E  A^V***-  A®  ,  E£*VP°‘>vpa  “B01  M 

vp  n  hsl  nlm  nu*  vp  n  hmc  mu*  nu*  (3-6) 

Furthermore,  since  we  work  with  unitary  irreducible  representations,  as 
shown  by  the  structure  of  the  operators,  aquation  (3-3),  tha  coefficients 
a  have  the  property 
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Now  wm  have  to  decoapoae  the  incident  field  into  parts  belonging  to  tha 
rows  of  irreducible  rcpraaentatlons  of  tha  symmetry  group  of  the 
cluster.  This  is  easily  done  for  tha  projection  operators,  aquation  (3- 
3),  have  the  completeness  property15'16 


SISDONI 


Thus  we  can  write 
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(3-8) 


(3-9) 
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and  an  analgous  equation  for  iB  ,  where 


Jvi»f  S  Os>Mkr>W:> 


and 


"§e2A(kr)?u. 

L«"5  ;ir^xjL(kr)XLM-(,1)- 


This  completes  the  syranetrization  procedure  because  equation  (2-4)  need  not  be 
syi.inietrized  for  reasons  that  will  become  clear  in  the  next  section. 

4.  EQUATIONS  FOR  THE  COEFFICIENTS 

We  are  now  able  to  write  the  equations  for  the 
coefficients  4^  and  of  the  scattered  wave,  equation  (3-5),  by 

imposing  to  E  and  d  the  appropriate  boundary  conditions  at  tne  surface 
of  each  one  of  the  spheres.  To  this  end  we  rewrite  equations  (3-5)  and 
(3-9)  in  terms  of  vector  harmonics  centered  at  a  single  site,  say  t^, 
through  the  use  of  the  appropriate  addition  theorems.  7 Indeed,  near 
the  surface  of  the  a-th  sphere  belonging  to  the  o-th  set  we  have 
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and  analogous  expressions  for  iiP8'  and  .  In  equation  (4-1)  we 

define  ~n  ~n 


L'M*LM 


^■^ag)£C(l,i-'  ,l'  ;-u,w+u)4it£  1l  l  *1  (l'  , m’ +u;l,m+u) 


X  hx(kRa6)VtM<-M^a6)C(l.t,L;-u,M+u) 


(4— 3a ) 


^wiru.  *  l;-y.ir*w)4w2:i‘'  *t.A*i 


x  Ix(*.-i,M-+M;L,efw)hx(Wa6)y*M._M(Ra6)C(l,t,L;.u,M+u) 


(4-3b) 


with  J^°e  -B  -?a  .  while  in  equation  (.4-2)  fywui  and  ^l-m-lm  are  identical 

to  and  respectively  but  for  the  substitution  of  jA  to  hA 

and  R6  *0 .  The  quantities 


VL* 


M' ;LM)  «  . /y* 


M* 


are  the  well-known  Gaunt  integrals.1^ 

How,  taking  the  dot  product  of  equations  (4-1),  (4-2)  and  (2-4)  in  turn 

With  ~cl  Y2m  (-£a}  and  x  <£«>  we  set  tha  radial  and  the 

tangential  components  of  the  field  at  the  surface  of  the  a-th  sprtere. 

Imposition  of  the  boundary  conditions  and  integration  over  the  angles 

yields,  for  each  t,m  six  equations  among  which  the  coefficients  of  the 

internal  field,  Cy^  and  D®  ^  ,  are  easily  eliminated.  This 

circumstance,  on  one  hand n  clarifies  the  inessentiality  of  the 

symmetrization  of  the  internal  field,  and  on  the  other  hand  yields,  for 

each  v,  p,  rQ  a  system  involving  only  the  "'s  and  s  's  as  unknowns: 


eee!ee(4  4 « r*flr+*a8 

T  8GT  L  <  M  m€n  v  Ll  **"  ^Sl.  <■  J  imLM  /  MUW  ml  X$GiK  tmLa  WJ*  iwrt  I 
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where  we  define 


\ikrj3T<rSi>  ■  ,taf,rJ,lkr,» 


(4-5a) 


r  »b 
a  a 


-  •>Fdjr.w> 

\ikr.>3^rJt>  ■  ’Vxb'MK » 


r  »b 
a  a 


(4-5b) 


Since  the  spheres  wichin  a  sec  are  idencicai  to  each  other,  and  do 
not  actually  depend  on  a  but  only  on  the  set  igdex,  Therefore, *from 

now  on  these  quantities  will  be  indicated  by  and  S°  ,  respectively. 
Equation  (4-4)  can  be  put  in  a  more  symmetrica*  form  through  the  use  of  - 
equation  (3-6).  Indeed,  multiplication  of  equation  (4-4a)  by  (a^) *  , 

and  of  equation  (4-4b)  by  (•>)!&„)*  *  and  summation  over  the  a's  belonging 
to  the  o-th  set  and  over  m  yield  for  each  v,  p,  t  the  equations 


Y  ™  L  n  n*  I’  r*’  ”*  n 

e(5  (e))sVP<,s^v  (e»m)A  p  *  -q  p 

r  VvL»J  <•’»•  /■»»•«  ’>•  ”* 

In  equation  (4-7)  we  put,  for  the  sake  of  simplicity, 

r  =  (a,*,n),  s  =  (t.l.n),  r’  =  ( a,i,n ’),  s'  =  (t.l.n-  ) 


(4-7*) 


(4-7b> 


and  define 
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'*  am  SM  nlm  lmU*  "LM 

xv(»,.)  -  E  L(»T)*Kae  6vpB 

«  am  SM 


(4-8) 
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with,  an  obvious  meaning  of  the. Decameters  e*..au...The  quantities  .  r$' 
and  cj!,("rf)are  Identical  to  1  a^(m)  and  cj^(#,m)  respectively,  except  for 
the  mutual  exchange  of  the  ars  with  the  b's. 

Moreover  - 1 - - - - 

jj  Err,  (ii  Ini’*  ,*  /•  1 1.*’*') 

—  ua  a.  w  nui  L  iwr  tmur  um  tour  1  \  *to  / 

Q*  -ZEE  “  (k)fcvp  4a  ♦ndvp<“  1  /bvpaV 

nr  lm  am  *  nU*  '  l  LMM*  Lm'Mr  lmw  tmurj  \  "ImJ 
We  remark  that  in  equations  (4-7)  and  (4-d)  the  superscript  p  on 
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ii  and  K  is  missing.  As  will  be  shown  later,  these  quantities  are  ' 
actually  independent  of  the  row  index  so  that  p  has  been  dropped. 

5.  THE  CROSS  SECTIONS 

Once  the  coefficients  %  and  are  known,  all  the 

scattering  properties  of  the  cluster  can  easily  be  calculated.  For  this 
purpose,  symmetry  is  not  essential  and  we  can  revert  to  the  unsyometrized 
expansion  of  the  scattered  field  through  equation  (3-6).  It  is 
convenient  to  express  the  scattered  field  in  terms  of  vector  harmonics 
centered  at  a  single  point  through  the  addition  theorem  already  used  in 
the  preceding  section. If  E  is  the  vector  position  of  the  chosen 
point  and  r  •  r  -  R 

'WQ  **  - - - - - - 
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l-mv*  an<*  LMLM'  are  given  by  equation  (4-3)  but  for  the  substitution 
of  to  .  Equation  (5-1)  is  valid  provided  that  r#>  ■  |  , 

i.e.,  in  the  region  outside  a  sphere  centered  at  R  ann  including  the 
whole  cluster.  Therefore  choosing  K  ■  0  and  thus  letting  the  center 
of  the  expansion  (5-1)  coincide  wit?j  the  center  of  symmetry  of  the 
duster,  the  .radius  ^of  the  above  sphere  is  minimized.  Anyway  the 
coefficients  A  and  B,  unlike  those  of  the  field  scattered  by  a  single 
sphere,  depend  on  the  direction  of  the  incident  wavevector,  k  .  As  a 
consequence  all  the  quantities -of  interest  depend  both  on  k  and  on  the 
scattered  wavevector,  k  “  kr,  except,  of  course,  the ~ scattering, 
abeorbtion  and  total  cross-sections  which  depend  only  on  k.  A 
straightforward  calculation  shows,  in  fact,  that  ^ 
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Finally,  we  notice  that  the  cross  sections  depend  on  the  polarisation  of 
the  incident  wave,  n,  as  explicitly  indicated  in  equation  (5-3). 

6.  DISCUSSION 

In  the  preceding  sections  we  used  group  theory  to  put  the  system  for 
the  coefficients  of  the  scattered  wave  in  factorized  form.  Let  us  now 
rewrite  equation  (4-7)  inatrixwise  as 


s"1  + 

K>,e)  4vp  Pvp 

5-1  ♦  »v(s>  avp  Qvp 

in  self  explanatory  notation.  Clearly  the  dependence  on  p  of  the  right- 
hand  side  of  equation  (6-1)  forces  us  to  solve  all  the  system  arising 
froa  the  factorization  procedure.  Nevertheless,  the  matrix  on  the  left- 
hand  side  of  equation  (6-1)  does  not  depend  on  p,  just  as  in  the  case  of 
a  secular  problem.  The  p-independence  of  R  and  s  ,  indeed,  follows  from 
their  very  definition,  equation  (4-5).  These  quantities  occur  even  in 
the  theory  of  electromagnetic  scattering  from  a  single  sphere^  and  are 
the  electromagnetic  analog:  of  the  elements  of  the  transition-matrix  in 
quantum  scattering  theory.  Rv  and  ,  in  turn  are  the  symmetrized 
counterparts  of  and  ,  equation  (4-3),  which  are  the  matrix  elements  in 
the  site  and  angular  momentum  representation  of  G,  the  dyadic  Green's 
function  for  free  space  propagation  of  spherical  vector  waves.  The  p- 
independence  of  R  and  *  is  then  a  direct  consequence  of  the  invariance  of 
G  under  the  symmetry  operations.  Therefore,  the  matrix  on  the  left  of 
equation  (6-1),  which  can  be  interpreted  as  the  inverse  of  the  transition 
matrix  for  the  whole  cluster,  2  turns  out  to  be  independent  of  p.  This 
circumstance  greatly  reduces  the  computational  work  in  the  case  of 
multidimensional  irreducible  representations.  Anyway  the  computational 
work  needad  to  aolva  our  scat taring  problem  depends  on  the  number  of 
spheres  in  the  cluster  and  on  the  number  of  L-values  included  in  the 
expansion  of  the  scattered  field,  equation  (2-3)  or  (3-5).  For  a  duster 
of  N  spheres  and  including  multlpolea  up  to  the  .order  Ljj,  the  dimension 
of  the  unsymmetrised  system  is  d  ■  2N(Lu  +  I)2  -  2N,  while  if  group 
theory  is  used,  the  order  of  the  factorized  systems  depends  on  the  group 
of  the  cluster.  As  an  example  1st  us  consider  a  cluster  of  5  spheres 
with  point  group  Tj  (the  SO  4,  CH4  molecules  have  just  this 
structure).  Table  1  snows  the  order  of  the  symmetrized  systems  and  that 
of  the  unaymmstrlzed  one  for  values  of  L  up  to  4.  Anyway,  from  the  work 
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of  hie^  and  of  Debye^we  know  that,  for  a  sphere  of  radius,  b,  *  1 
is  quite  sufficient  for  a  good  description  of  the  scattered  field, 
provided  that  Kb  «  1.  Wow,  since  our  cluster  is  meant  to  modelllze  a 
molecule,  the  radius  of  the  spheres  can^hardiy  exceed  3  A  so  that  in  the 
visible  and  Infrared  range  Kb  <  10  for  any  a.  Therefore,  even 
allowing  for  the  effect  of  multiple-scatterings,  •  3  should  be  quite 
sufficient  to  achieve  fairly  converged  values  of  the"  scattered  field. 

Table  1 


I’M 

i 

2 

3 

HI 

n 

i 

2 

6 

10 

i 

2 

6 

10 

H 

2 

a 

12 

20 

H 

4 

10 

19 

30 

m 

4 

10 

19 

30 

H 

30 

- 

150 

240 

Dimension  of  the  symmetrised  and  unsycuaetrized  systems  for  up 
to  4  for  a  cluster  of  5  spheres  with  a  point  group  Td.  The  entry  U  means 
unsynunetrized  while  the  other  entries  indicate  the  irreducible 
representations  in  the  notation  of  Haiaernesh. 

7.  COMPUTER  CODES 

The  method  described  has  been  implemented  through  computer 
algorithms  whose  practical  application  we  briefly  describe.  The  first 
code  has  a  main  program  named  SYMULT  and  subroutines  PRMUTE,  PRJECT, 
FUNCTION  DL,  PKJPW,  0RTH0G.  The  second  main  program,  named  MUSCA,  has 
subroutines  RESYCO,  REFAPW,  GHIlT,  GJMT,  GHGK,  G1GL,  SHHME,  SKMilE,  SPQ, 
MULAM,  CROSSE,  CMES,  ZOTS,  WLM,  RSINV,  AUXIL,  FUNCTION  DCLEB,  RBF,  RNF , 
CBF,  POLAR,  SPHAR,  DC ML IN,  COMPLEX  FUNCTION  DCDOT. 

SYMULT  only  reads  the  input  data  and  calls  the  subroutines, 
which  actually  parfora  the  aatheaatical  computations,  and  finally  prints 
the  output,  which  consists  of  the  following  quantities: 

i.  Names  of  the  operations  for  the  group,  ,  with  the 
corresponding  pereaetere:  u  ,  angle  of  rotation,  X  ,  p,  v  direction 
coalnes;  end  the  permutations  induced  on  the  sites  of  a  given  set  by  the 
operations  of  .  For  the  operations  of  group  and  associated  operator 
ccc  we  refer  to  eqs*  (3-2a)  (3-2b)  end  following  explanations. 

2*  Coefficients  of  the  independent  combinations  of  magnetic  and 
electric  nultlpoles  centered  at  the  sites  of  the  given  set,  such 
quantities  ere  computed  from  eqs.  (3-4a),  (3-4b). 
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3.  The  projeccions  of  the  magnetic  and  electric  muicipoies 
included  in  the  expansion  of  the  incident  plane  wave,  as  are  shown  in  eqs 
(3-9). 

These  quantities  are  written  on  a  file  and  are  the  input  data 
for  MUSCA,  whose  output  consists  of  the  following  quantities: 

1.  The  total  cross  section,  computed  from  eq  (,3-5c) 

2.  The  scattering  cross  section,  from  eq  (3-5a) 

3.  The  absorption  cross  section,  from  eq  (3-5b) 

4.  The  coefficient^  of  the  scattered  wave  for  each 
representation,  the  Anfi£  and  BnfiE  of  eq.  (3-6),  which  are  the  unknown 
quantities  in  our  solution  and . thus,  with  their  knowledge  we  are  able  to 
fully  describe  the  properties  of  the  cluster,  as  stated  in  section  5. 

Following  are  some  (lata  obtained  from  these  programs  (outputs  of 
SYMULT  and  MUSCA)  for  the  case  of  the  point  group  Td  (SO^  )  and 
comparison  is  made  with  experimental  data. 
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